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Abstract 

Topological defects are common in many everyday systems. In general, they appear if a 
symmetry is broken at a rapid phase transition. In this article, I explain why it is believed that 
they should have also been formed in the early universe and how that would have happened. 
If topological defects are found, this will provide a way to study observationally the first 
fractions of a second after the Big Bang, but their apparent absence can also tell us many 
things about the early universe. 
DAMTP-2003-71 

1 Introduction 

Imagine sitting at a round dinner table, with a plate in front of each guest and a spoon at half 
way between each two plates. This setting is symmetric between left and right, and you could 
therefore choose to use either the spoon on your left or the one on your right. However, as soon 
as your neighbours make their choice, this symmetry disappears, and you have to conform to that 
choice. This is a classic example of spontaneous symmetry breaking (SSB) (see Ref. PQ). 

However, it may happen that the people on your right choose to use the spoons on their left 
and vice versa, leaving you with no spoon at all. This is an example of a topological defect. They 
occur generally in systems with SSB, because that requires a choice between several identical 
possibilities, and the choice is often different in different parts of the system. 

Spontaneous symmetry breaking is a very common phenomenon in many physical systems. 
In the Standard Model of particle physics, it is in a sense responsible for the masses of the 
elementary particles, although the associated symmetry is of a more abstract nature, as it is a 
local gauge symmetry rather than a global symmetry. The Standard Model actually happens to be 
a special case with no topological defects, but this, nevertheless, suggests that its extensions would 
generally predict the existence of topological defects. Depending on the details, these defects could 
be domain walls, cosmic strings or magnetic monopoles. 

Both in the early universe and in our dinner table example, the symmetry was initially un- 
broken. It got broken when the first guest picked up a spoon, or when the temperature of the 
universe decreased below a certain critical temperature. At the dinner table, it is easy to see that 
if the table is large enough, topological defects are bound to form but the process of defect 
formation was more complicated in the early universe. The purpose of this article is to explain it. 

We can also ask if these defects could have survived until today and if and how they could be 
observed in that case. So far, we have no observational evidence for their existence. This alone 
can tell us many things about the early universe. On the other hand, defects have not been ruled 
out, and if they are found, that would have a huge impact on our understanding of the universe. 

In order to draw any conclusions from the absence or possibly the existence of defects, it 
is crucial that we understand the process of defect formation properly. Fortunately, we do not 
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Natural units 


SI units 


Energy 


1 GeV 


1.60 x 1CT 1U J 


Temperature 


1 GeV 


1.16 x 10 13 K 


Mass 


1 GeV 


1.78 x 1CT 27 kg 


Distance 


1 GeV" 1 


1.97 x 1CT 16 m 


Time 


1 GeV" 1 


6.65 x 1(T 25 s 



Table 1: Conversion between natural and SI units. 



have to rely solely on theoretical calculations, because topological defects are also formed at 
phase transitions in certain condensed matter systems such as superfluids and superconductors. 
This phenomenon is theoretically very similar to its cosmological counterpart, and we can use 
this analogy to do 'cosmological experiments' On a more fundamental level, these same 
experiments can be used to test our understanding of non-equilibrium dynamics of quantum field 
theories, which will be very important for particle physics. 

In this article, I review the theory of defect formation at phase transitions, emphasizing those 
aspects that are relevant for the early universe. I will also discuss the importance of defects to 
cosmology, but readers who want to know more about that are advised to read Refs. ^JOOH)]. 
More discussion on defect formation at phase transitions can be found in Refs. 0121 ED- 

Throughout the paper, I will use the natural units c = h = Ub = 1. This means that everything 
is expressed in units of GeV (giga electron volt). Tabled shows how they can be converted into 
SI units. 

2 Scalar fields and global symmetries 

Apart from gravity, all physics at microscopic scales is described by quantum field theories (see 
Ref. 10 ). For the moment, we will restrict our discussion to scalar fields, which is the simplest 
type of quantum fields. A scalar field is Lorentz invariant, which means that its value does not 
depend on the reference frame. It simply has some value at each point in spacetime, and this value 
is independent of the orientation and velocity of the observer. These values can be real numbers, 
but we can also think of complex scalar fields, or vector or matrix valued scalar fields. 

Because the field is Lorentz invariant, there is a large amount of freedom in constructing scalar 
field theories. Most of this freedom is in the choice of the potential V(cf>), which gives the energy 
density of a constant field with value (f>. For a real scalar field </>, the potential leads to the classical 
equation of motion 

<^-v 2 = -y'(0), (i) 

where the double dot indicates a second time derivative, and the prime indicates a derivative with 
respect to (j>. Apart from the gradient term V 2 4>, Eq. is identical to the Newton equation for 
a particle in the same potential V. Therefore, we can understand the behaviour of the scalar field 
by thinking of a ball moving in a one-dimensional landscape. 

In the vacuum state, the field is constant in space and time, and its value corresponds to the 
global minimum of V(<p). The value of <f> in the vacuum is known as the vacuum expectation value 
(vev). 

Let us now imagine a potential that has reflection symmetry around 4> = 0. Typically, this 
would be a polynomial in </> 2 , for instance 

V{<j } ) = \x^-v 2 ) 2 . (2) 

The value of the potential does not change if we flip the sign of <f). If v 2 is negative, the vacuum 
state is at = 0. However, if v 2 > 0, this point becomes a local maximum, and two minima 
appear at 4* = i". The field is now in the same situation as a guest at the dinner discussed at 
the beginning of this article. It will have to choose one of the two minima, and this breaks the 
reflection symmetry spontaneously. 
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Figure 1: A U(l) symmetric potential for a complex scalar field. The ball depicts the value of 
the field in a vacuum state at the bottom of the potential. There is a circle of possible vacuum 
states, characterized by the phase angle 6 — arg cf>. 

If we have, instead of a real field, a complex one, we can think of a ball on a two-dimensional 
surface, which corresponds to the complex plane. We will be mostly interested in the case in 
which the potential depends only on the absolute value of cj>, not on its phase angle. Then, this 
two-dimensional picture is symmetric around the origin. A potential like this is shown in Fig. ^ 
A symmetry with respect to rotations around one axis, such as this, is known as U(l). 

We can easily think of scalar fields with more components, although it is then more difficult to 
visualize the motion. In that case, the mechanical analogue is a ball moving in a higher-dimensional 
landscape, which we will call the internal space. 

Whenever a symmetry is broken, the potential has several minima, which are related to each 
other by the symmetry, and each of them corresponds to a possible vacuum state. If the symmetry 
is continuous, these minima form a continuous valley with a perfectly flat bottom. If we give the 
ball a push in this direction, it will roll around the set of possible minima, no matter how weak 
the push was. This has the consequence that the quantum theory has states with arbitrarily low 
energies, or in other words, massless particles. These particles are known as Goldstone bosons. 

Our previous dinner table example is not very helpful for visualizing the breaking of a U(l) 
symmetry, because it only dealt with a choice from two possibilities. In the case of U(l), we have 
a continuous set of possible vacua. Perhaps the closest analogue in everyday life is the time of 
the day. If we did not interact with other people, it would not matter to which time we set our 
watches. At any time of the day, we could choose any possible value from the whole range of 
possibilities from 00:00 to 23:59. Because 24:00 is the same as 00:00, we can think of a circle of 
possibilities. This can be interpreted as a U(l) symmetry. 

However, our interaction with other people breaks this symmetry. The whole society has to 
set their watches to the same time in order to work properly, and this time plays the role of the 
vacuum expectation value. In the UK it was chosen that the clocks are set to 12 o'clock when the 
sun is at its highest above Greenwich, but in principle this choice was arbitrary. 
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3 Gauge fields and local gauge symmetries 



The symmetries we have discussed so far have been global symmetries: The system is only sym- 
metric under rigid rotations, where the field is rotated by the same amount everywhere. In our 
time analogy, this corresponds to changing the time by the same amount in the whole country. 
This is, of course, done every spring and autumn, when the UK changes from Greenwich Mean 
Time (GMT) to British Summer Time (BST) and back, and because of this symmetry, it does not 
cause any problems. On the other hand, if you forget to move your clock, you will most probably 
be missing appointments and encountering other problems. 

The Standard Model of particle physics actually possesses an higher symmetry that allows 
rotations of the fields at different points by different amount without changing their physical 
content. A rotation like this is known as a local gauge transformation, and the corresponding 
symmetry as local gauge symmetry or gauge invariance. It is made possible by the gauge fields, 
which also carry the electromagnetic, weak and strong interactions. In this way, this symmetry 
dictates many properties of these interactions. 

We can use the time analogy to illustrate gauge fields. As we noted, you cannot use a separate 
time convention from the rest of the country without having problems, but worldwide, different 
countries do use different time conventions. This is made possible by an extra structure, namely 
the time zones. We know that when we travel from the UK to, say, New York, we have to move 
our watches five hours backwards, and this is analogous to what the gauge field tells us in the case 
of a gauge field theory. If we want to compare the scalar field to its value at some other point, we 
have to first rotate it by the amount indicated by the gauge field. 

Every March, when the UK moves clocks one hour forward and changes to BST, we carry out 
something similar to a local gauge transformation. Many countries do not move their clocks at 
the same time, but that does not cause problems as long as travellers are aware what the current 
time difference is on any given day. Analogously, if we carry out a local gauge transformation in a 
gauge field theory, the gauge field is changed in order to compensate for the rotation of the scalar 
field. 

The simplest example of a gauge field theory is the Abelian Higgs model. It consists of a 
complex scalar field (j> with a U(l) symmetry and a gauge field A, which is a real vector-valued 
field. 1 Because we cannot directly compare field values at two different points, we have to replace 
gradients by covariant derivatives 

Vcj> -> = + ieA(f>, (3) 

where e is known as the gauge coupling. (In non-Abelian theories, which we will discuss later, the 
gauge coupling is usually denoted by g.) Physically, this replacement gives the scalar particles an 
electric charge e, with the gauge field playing the role of the vector potential. The electric field is 
given by E = —dA/dt, and the magnetic field by B = V x A. 

In the broken phase, where (f> has a non-zero value, the system becomes essentially a super- 
conductor. Any field configuration with a constant and non-zero <f> and B would have an infinite 
energy density. This is known as the Meissner effect and is a characteristic property of supercon- 
ductors (see Ref. JI]). Real superconductors can be described with a similar model, in which the 
role of 4> is played by Cooper pairs that consist of two electrons. In particle physics and cosmology, 
we are would think of the Abelian Higgs model as a relativistic theory and as a fundamental 
field. Then, the Meissner effect corresponds to a non-zero photon mass and is also known as the 
Higgs mechanism. In contrast to the global case, there are no massless Goldstone bosons, because 
the rotation around the vacuum manifold can always be compensated by the gauge field. It is 
therefore often said that the gauge field "eats" the Goldstone boson. 

In terms of our time zone analogy, the presence of a non-zero magnetic field would correspond 
to a hypothetical situation, in which the amount by which travellers need to move their clocks 
when they cross a border is not determined globally by time zones, but is specified for each border 

1 Actually, the gauge field is a four- vector Au, but we ignore the time component here. This can be done 
consistently without changing the physical content of the theory, and is known as the temporal gauge. 
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separately. Imagine that when you travel from the UK to France, you are told to move you clock 
one hour forward. When you travel further to Spain, you are again told to move your clock forward 
by one hour. Finally, when you travel from Spain to the UK, you are told to move your clock one 
hour backward. As a net effect, when you return home, your watch would show one hour more 
than the clocks you have in your home. In reality, this does not happen because each country 
follows a specific time zone. To the extent that the choice of the time zone corresponds to SSB, 
we can understand this as an example of the Meissner effect. 

4 Particle physics 

In particle physics, we encounter more complicated gauge symmetries than U(l). Weak and strong 
interactions are described by groups knowns as SU(2) and SU(3), correspondingly. These groups 
are called non- Abelian, because in these cases the gauge field is matrix valued. The product of 
two matrices depends on the order, i.e., AB ^ BA, and in algebra, an operation with this propcrt 
is called non-Abelian or non-commutative. Nevertheless, the basic principles are the same as with 
the Abelian U(l) group. 

The Standard Model of particle physics is based on SU(3), SU(2) and U(l) groups, and gives 
an extremely accurate description of strong, weak and electromagnetic interactions in terms of 
gauge fields, quarks, leptons and a scalar field known as the Higgs field. 

The Higgs field is a complex two-component scalar field, which is sensitive to both U(l) and 
SU(2) groups. In the vacuum, the Higgs field has a non-zero value (v w 246 GeV), which breaks 
these groups spontaneously, but a residual U(l) subgroup survives. This U(l) group gives rise to 
electrodynamics in the same way as in the Abelian Higgs model. As a consequence, the Higgs 
mechanism makes the W and Z bosons massive, but leaves the photon massless. 

The Standard Model has been highly successful experimentally, and in fact, there are still no 
experimental results that disagree with it. Nevertheless, it is in a sense unnecessarily complicated, 
because it has three separate gauge groups, each with its own coupling constant, and the way the 
quarks and leptons transform under gauge transformations does not seem very natural. In fact, if 
one calculates how the theory behaves at high energies, one finds that at around 10 16 GeV the 
strong, weak and electromagnetic interactions become equally strong. Furthermore, if one thinks 
of these gauge groups as subgroups of a larger gauge group, it turns out that it is also possible 
to explain the gauge transformation properties of quarks and leptons. There are several possible 
ways of doing this, and they are known as grand unified theories (GUT). The simplest example is 
based on a group known as SU(5). 

In GUTs, the grand unified gauge group is split into the SU(3), SU(2) and U(l) subgroups 
by SSB. In the SU(5) theory, that is achieved by postulating a new, 24-component scalar field <&, 
which has a non-zero vev (of the order 10 16 GeV). Unfortunately these energies are far beyond the 
reach of any experiments. There are other GUTs in which the U(l) group appears only at lower 
energies, such as the SO(10) theory with the symmetry breaking pattern 

SO(10) -> SU(4) x SU(2) x SU(2) 

-► SU(3) x SU(2) x U(l), (4) 

where the energy scale of the latter symmetry breaking could be of the order 10 12 GeV |13| . 

5 Topological defects 

So far, we have been mainly discussing the vacuum states, but it is often also possible to find other 
time-independent states. Topological defects are one class of these, and they arise in models with 
spontaneously broken symmetries, if the choice of vacuum is different in different directions and 
these choices do not match completely. 

We already encountered one very simple example of a topological defect in the introduction. 
There, the defect arose because the symmetry breaking involved a choice from a discrete set of 
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possibilities. Defects like these are called domain walls, because in three dimensions, such a defect 
would form a surface separating regions of different vacua. For a real scalar field with the potential 
in Eq. @ the corresponding solution is the "kink" , 

<fi(x, y, z) — v tanh \f\vz. (5) 

This solution has a domain wall at z = 0, where the field vanishes and the symmetry is restored. 
On different sides of the wall, the fields approach different vacua asymptotically. The kink is stable, 
because it is not possible to deform the field into the vacuum configuration without changing its 
value everywhere in half of the space simultaneously. However, if a kink encounters an anti-kink, 
they annihilate. 

In the case of a U(I) symmetry, the set of possible vacua, i.e. the vacuum manifold M, is a 
circle (see Fig. I). Let us now consider a closed curve, which we will denote by C, in space. At 
each point on the curve C, the field <p has some value, which we can think of the position of a ball 
in the potential in Fig. I. Let us now move along the curve C and ask how the ball moves in the 
potential. 

In vacuum, <p has the same value everywhere, and the ball would therefore simply stay still in 
the same place, but this is not always true. In any case, when we have moved around the whole 
curve C and come back to our starting point, the ball must also have returned to its original place, 
but along the way, it may have travelled around the circle of minima. The number of times it 
moved around the circle is known as the winding number. 

Let us now assume that this winding number is not zero. If we now deform the curve C 
continuously, the path of the ball in the potential must also change continuously. If the ball 
stays on the vacuum manifold, the winding number cannot change, because that would require 
a discontinuous change. On the other hand, if wc shrink the curve C to a point, the path can 
consist of one point only and therefore must have winding number zero. The only way to avoid 
a contradiction is that the path of the ball has left the vacuum manifold at one point, but this 
means that somewhere inside the curve C, there must have been a point at which (j) was not in 
vacuum. This point is known as a vortex. 

In circular cylindrical coordinates, a typical vortex has the rotation-invariant form 

4>{r,<p,z) = ve i * n ™f(r), (6) 

where riw is the winding number and /(r) is a function that vanishes at r = and approaches 1 
at infinity. In three dimensions, the vortices are one-dimensional vortex lines (see Fig. |2J, and in 
the cosmological context, they are also known as cosmic strings. 

We have seen that the time zones correspond to a spontaneous breakdown of a U(l) symmetry 
and there are also vortices in this analogy, namely the North and South Poles. Around each pole, 
the time zone changes by a full 24 hours around any closed path that encircles a pole. The time 
zone is not well defined at the pole itself, and the symmetry is therefore restored there. 

For a three-component real scalar field, the vacuum manifold is a sphere. In that case, it is 
possible that the field values on a closed surface S in space wrap around the sphere, analogously 
to how the field values wrapped around the circle on the curve C if the winding number was 
non-zero. Then the same argument shows that there must be a point inside the surface at which 
the field leaves the vacuum manifold. This is a monopole. A typical monopole has the spherically 
symmetric "hedgehog" form, 

$ Q (x)=^/(|x|)^, (7) 

where v is a constant. The function f(r) is again a function that vanishes at r = and approaches 
unity at infinity, which means that the field is on the vacuum manifold far away from the origin. 
In $ a the index a indicates the direction in the internal space, whereas in x a it indicates the 
direction in coordinate space. In this sense, <E> and x are parallel. 

With a four-component scalar field, there is a further possibility that the field configuration in 
the whole three-dimensional coordinate space wraps around the vacuum manifold. That is known 








as a texture ^3]. It is easy to see that a texture is not a stable object, because the energy of any 
scalar field configuration is generally of the form 

E = Jd 3 x [(V0) 2 + Vffl] . (8) 

In a texture, the field is everywhere on the vacuum manifold, and the latter term vanishes. The 
energy of a texture is completely due to the spatial variation of the field (f>. If we double the size 
of the field configuration, the first term doubles, and therefore a texture would shrink to a point 
and disappear in order to minimize its energy. 

In fact, only domain walls are truly localized objects in global theories, because the energies 
of vortices and monopoles diverge logarithmically and linearly with the system size, respectively. 
This is again due to the gradient term in Eq. © . The only way to have a configuration with a finite 
energy is to have an equal number of vortices and antivortices, or monopoles and antimonopoles. 
Even then, a vortex and an antivortex would have a logarithmic interaction, which binds them in 
a pair. For a monopole-antimonopole pair, the interaction is linear, and therefore they would be 
confined just like quarks. 

Topological defects also exist in gauge field theories, but their properties are somewhat different. 
Far from the defect, where the scalar field is on the vacuum manifold, the gauge field can cancel 
the gradient contribution to the energy in Eq. (jHJ), because the gradient V<^> gets replaced by the 
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covariant derivative T)(j> (see Eq. ©). In the case of the Abelian Higgs model, this happens when 



where 9 = arg <f> is the phase angle of the scalar field. Thus, the (cross-sectional) energy of the 
vortex is finite, but it still does not vanish completely, because the energy density is non-zero near 
the vortex core, where the field is away from the vacuum manifold. In any case, this means that 
the energy of a defect is localized and finite in gauge field theories. Their interactions are also 
much weaker than those of global defects: Interactions of vortices are exponentially suppressed at 
long distances, and monopoles have a Coulomb type interaction. 

It also follows from Eq. that a gauged vortex carries magnetic flux Q2J. The phase angle 
changes by an integer multiple of 27r around a vortex, but this means that the contour integral 
of the gauge field around the vortex must be a multiple of 2ir j e. This contour integral is nothing 
but the magnetic flux through the contour, and therefore the flux is proportional to the winding 
number 



In other words, each vortex carries an integer multiple of the flux quantum $o = 2ir/e. 

This can be seen by placing a superconductor in an external magnetic field. As observed earlier, 
the superconductor repels the magnetic field, but in many cases, if the field is strong enough, it 
will penetrate the superconductor and form Abrikosov flux tubes, which are exactly what we have 
called gauged vortices. Superconductors for which this happens are said to be of Type II. In 
Type I superconductors, the magnetic field forms one thick vortex that carries the whole magnetic 
flux. If A < e 2 , the Abelian Higgs model behaves like a Type I superconductor. Vortices with 
multiple windings are stable and two vortices of equal sign attract each other. In the opposite 
case with A > e 2 , two equal-sign vortices repel each other, and the theory resembles a Type II 
superconductor. 

The simplest example of a gauge theory with monopoles is the Georgi-Glashow model. It 
consists of a real three-component scalar field & a , where a G {1,2,3}, and a three-component 
gauge field A a . The gauge field makes the theory invariant under local (i.e. position-dependent) 
rotations of <E> a in the three-dimensional internal space. This symmetry is spontaneously broken 
by a vacuum expectation value $ a $ a = v 2 . However, it is still possible to do rotations with <E> a 
as an axis without changing it, and therefore there is a residual U(l) symmetry. We have seen 
before that a local U(l) symmetry gives rise to the electromagnetic field, and this happens in this 
case, too. 

In normal electrodynamics, the field lines of the magnetic field cannot end, which means that 
there cannot be magnetic charges. The same is true in the broken phase of the Georgi-Glashow 
model as long as $ a ^ 0. However, $ a vanishes at the centre of the monopole solution in Eq. Q, 
and the corresponding solution in the Georgi-Glashow model does indeed have a non-zero magnetic 
charge. This solution, which is known as the 't Hooft-Polyakov monopole |lfil I17j . is therefore a 
magnetic monopole. The total energy of a 't Hooft-Polyakov monopole is finite, but it has a long- 
range 1/r contribution. Because the monopoles have magnetic charges, it is easy to understand 
this as the magnetic analogue of the Coulomb interaction. 

Similar 't Hooft-Polyakov monopoles exist in essentially all theories with a residual U(l) sub- 
group. In particular, this includes all GUTs, and therefore they generally predict the existence 
of magnetic monopoles. The masses of these monopoles are determined by the energy scale at 
which the U(l) group appears. In the SU(5) GUT, the monopoles would be extremely heavy, 
m,M ~ 10 17 GeV, but for instance in SO(10), somewhat lighter monopoles with tum ~ 10 13 GeV 
would be possible (see Section 2J. 

6 Early universe 

It was observed by Hubble in 1929 that the light of distant galaxies is redshifted, and that the 
amount of redshift is proportional to the distance. The simplest explanation to this is that the 
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galaxies are moving away from us and the redshift is caused by the Doppler effect. Because the 
velocity is proportional to the distance from us, this does not imply that we occupy a special 
position at the "centre" of the universe, but that the universe is expanding with the same rate 
everywhere. 

In order to study the evolution of the whole universe, we assume the cosmological principle, 
which means that the universe looks the same at every point and in every direction, which seems to 
be a good approximation at very long distances. In technical terms, we assume that the universe 
is homogeneous and isotropic. 

This assumption, together with the observation that the universe is spatially flat JS| means 
that we can describe its expansion by two equations, the Friedmann equation and the conservation 
of energy, 



Here a(t) is the scale factor of the universe, the dot indicates a time derivative, p is energy density 
and p is pressure. The scale factor shows how the space stretches. The distance between two 
objects at rest at time t\ is longer by a factor of a(ti) / a(to) than it was at time to- As the Hubble 
parameter H — a/a is positive, the universe is expanding, and distant galaxies would indeed 
appear to be moving away from us. 

Conservation of energy l|12jl shows that as the universe expands, the energy density decreases. 
Consequently, the temperature of the universe has been higher in the past. If we know the equation 
of state, i.e. the pressure p as a function of p and T, we can use Eqs. and l|12() to follow the 
time evolution of the scale factor a arbitrarily far back in time. If we assume that the equation 
of state does not change dramatically, this analysis shows that the scale factor started from zero 
about 13.7 billion years ^Bj ago. We should not take this result too literally, because at that time, 
the energy density would have been infinite and we have no reason to believe that our assumptions 
were valid. 

In any case, we can quite confidently extrapolate the evolution of the scale factor almost all the 
way back to this initial singularity, which we can choose to be at t = 0. By observing the cosmic 
microwave background (CMB) radiation, we can actually see all the way to t ~ 380000 years |18| . 
when the universe became transparent. Although we cannot see beyond this "surface of last 
scattering", we can use theoretical calculations to find out what happened before that. 

Atomic and nuclear physics, and particle physics that has been tested and confirmed in ac- 
celerator experiments, can take us all the way back to i x 10~ 10 s, which corresponds to the 
temperature of 100 GeV. To reach ever earlier times, we have to make assumptions about physics 
that is out of the reach of any experiments. On the other hand, we can use the universe as the 
ultimate accelerator experiment, and calculate what kind of observational predictions follow from 
a given theory. This gives us a way of reaching much higher energies than particle accelerators 
ever can, albeit with the price that we cannot repeat the experiment. 

As we will see in Section symmetries that are broken at zero temperature are usually un- 
broken at high temperatures. If the temperature of the universe was initially high enough, the 
universe would have undergone a phase transition from the unbroken to the broken phase as it 
expanded and cooled down. We therefore expect that there were a number of symmetry breaking 
phase transitions in the early universe, and they may have led to defect formation. 

On the other hand, this hot Big Bang scenario cannot have been valid all the way to the 
initial singularity, and perhaps the most compelling reason is the horizon problem 19 . Different 
directions on the sky look very similar, and if we actually measure the temperature of the cosmic 
microwave background, it is uniform up to one part in 100 000 in different directions. Yet, this 
light was emitted only 380000 years after the initial singularity, during which time any information 
could only have travelled a distance of about 1 degree on the sky. Therefore, there is no physical 
mechanism that could have made different directions reach thermal equilibrium with each other. 
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Figure 3: The shape of the potential in Guth's original inflationary model. The field, represented 
by the gray ball, is initially in the local minimum, and the non-zero potential energy makes the 
universe expand with a constantly accelerating rate. Finally, the field tunnels through the potential 
barrier, and inflation ends. The non-zero potential energy is released during reheating, and the 
universe reaches thermal equilibrium. 

There are also other problems with the Big Bang scenario, such as the absence of magnetic 
monopoles, to which we will return later. 

In order to solve these problems, Alan Guth proposed the theory of inflation (see also 
Refs. Hnil2])j whose basic idea is that the universe was dominated by vacuum energy at some 
early stage. This can be achieved with a real scalar field, which is known as the inflaton and which 
has a potential with the shape shown in Fig. |3J The energy density p is mostly due to the kinetic 
energy 4> 2 /2 and the potential energy V(<j)). If the field is initially placed at the local minimum at 
<t> = 4>o, the time derivative 4> vanishes, and only the vacuum contribution V(<po) remains. Because 
(f> does not change, the energy density p is constant and Eq. (|ll|l implies that the scale factor a(t) 
grows exponentially. 

Inflation not only explains why the universe is so homogeneous, but also where all the structure 
we see in it today came from. The exponential expansion amplifies quantum fluctuations, and 
they act as seeds for structure formation and can be seen as fluctuations in the CMB. Recent 
measurements of these fluctuations |18| agree remarkably well with the predictions. 

In order to explain the observed homogeneity over the whole sky, inflation must have expanded 
the universe by at least a factor of around 10 20 , but after that, it must have come to an end in 
order to match with the ordinary Big Bang cosmology. Because energy is conserved, the energy 
density that was driving inflation, was released as radiation during this process. The dynamics of 
this reheating depends on the details of the inflationary model, but eventually the universe reached 
thermal equilibrium at the reheat temperature T r h. Its value is not known, but the asymmetry 
between matter and antimatter requires it to be at least 100 GeV. 

In Guth's original inflationary model did not describe reheating in a satisfactory way, but there 
are a vast number of other inflationary models that have been proposed to address this and other 
unattractive features of the original proposal. Most of them are based on the assumption of slow 
roll. In an expanding universe, the equation of motion |JJ| is modified to 

4> + 3H<j)-\7 2 cf> = -V'{(f>), (13) 
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Figure 4: a) A typical potential for slow roll inflation. During inflation (A), the field is in a 
part of the potential that is almost flat, and therefore it is moving very slowly. This gives rise 
to almost exponential expansion. Inflation ends when the potential becomes steeper (B). The 
universe reheats, and the inflation field ends up in the minimum of the potential (C). b) The 
potential for hybrid inflation. At any time, the waterfall field x is a t its minimum, and therefore 
the potential is effectively of the type shown in (a). At the end of inflation, the \ ~X symmetry 
is broken spontaneously at a phase transition. 

where H = a/a is the Hubble parameter. The extra term in this equation behaves as a friction 
term and slows down the motion of the field in an expanding universe. Just like an object falling 
in the atmosphere reaches a terminal velocity due to the air drag, the scalar field reaches a slow 
roll solution 

i V'{4>) 



3H 



(14) 



If the slope of the potential is small compared with the Hubble rate, the value of the poten- 
tial changes only very slowly. The energy density p is, again, almost constant, and we find an 
exponentially growing, inflationary, solution. 

In this setup, inflation can end smoothly if the potential starts to fall steeply at some value of 
</>, as shown in Fig.^J Perhaps the most elegant way of achieving that is by adding another scalar 
field x with the potential 



V{<f>,x) = \™ 2 tf + \g 2 ^x 2 + \* (x 2 



(15) 



The shape of this potential is shown in Fig. 0] and it is symmetric under x ~ y ~X- This model in 
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Figure 5: a) The effective potential in a continuous phase transition. Above the critical tem- 
perature T c , the origin <\> = is the only minimum, and the symmetry is restored. Two degerate 
symmetry-breaking minima appear at T = T c and the origin becomes unstable, b) The effective 
potential in a first-order phase transition. At very high temperatures, the origin <fi = 0, is the 
only minimum, but when the temperature goes down, two metastable symmetry-breaking minima 
appear. At the critical temperature T c , these minima are degenerate with the symmetric one. 
Below T c , the symmetric minimum becomes metastable and eventually decays through bubble 
nucleation. 

known as hybrid inflation |22j . 

We can see that if 4> has initially a large value and x is zero, the fields roll down slowly along 
the valley towards the origin and the universe inflates. However, when <j) reaches 4>q — ^^v/g, 
the "waterfall" field \ has to choose between the positive and negative branches, and this breaks 
the symmetry spontaneously. In either branch, the potential is much steeper, and therefore the 
slow roll assumption fails and inflation ends. The interactions between the inflaton and other 
fields reheat the universe as the inflaton rolls quickly down the hill to the minimum. 

There are many variants of this hybrid inflationary model. One attractive feature of the models 
in this class is that potentials of the form of Eq. i|15|) arise relatively naturally in many particle 
physics models. As we have seen, inflation is ended by a symmetry breaking phase transition in 
these models, and it is therefore important to understand what happens at this transition and, in 
particular, whether topological defects are formed. 

7 Thermal phase transitions 

At zero temperature, the equilibrium states of the system correspond to the minima of the po- 
tential. However, it is very well known from many everyday systems that a symmetry that is 
spontaneously broken at low temperatures can be restored at higher temperatures. One exam- 
ple of this is water: The crystal structure of ice breaks the rotation invariance, but when the 
temperature is above 0°C, ice melts and the rotation invariance is restored. 

In field theories, one can see this by defining the effective potential V c g (cf>), which is essentially 
the free energy of a system in which the spatial average of the scalar field </>(x) is constrained to 
be <j). Because the free energy is minimized in thermal systems, the equilibrium state corresponds 
to the minimum of the effective potential. 

The effective potential is typically calculated using perturbation theory [33], and the leading 
contribution comes from one- loop Feynman diagrams. In the scalar theory, the main effect is to 
give a temperature-dependent quadratic contribution to the potential, AV(<p) oc \T 2 (f> 2 |24l 125) . 
As shown in Fig. E^, this restores the symmetry when the temperature is above the critical 
temperature T c w v. 
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The minimum of the effective potential gives the mean value (0) , but the field fluctuates around 
this value. Therefore, even if the symmetry is restored and the mean value (</>) vanishes, the actual 
value of the field at any given point is typically non-zero. In fact, the field only vanishes in a set 
of measure zero at any given time. 

To leading approximation, the thermal fluctuations have a Gaussian distribution and therefore 
their properties are specified by the two-point function, which typically has an exponentially 
decaying form. The decay rate is characterized by the correlation length £. In practice, it means 
that fluctuations at two points separated by a distance greater than £ are independent. 

When the transition point is approached, the correlation length diverges with a critical expo- 
nent v 1 



which depends on the universality class of the system. At the same time, the dynamics of the 
system becomes slower, and this critical slowing down can be expressed in terms of the relaxation 
time t(T), which diverges with a (different) critical exponent /j,, 



For a real scalar field v w 0.630 and for a complex field v w 0.672 [27]. In most systems, 
/' 1- 

Unfortunately, it is often not possible to use the simple description in terms of the effective 
potential in gauge field theories. Thermal fluctuations can easily change the field configuration to 
any of its gauge transforms, because that does not cost any energy. This means that if we average 
over these fluctuations, the mean value (cj)) is always zero. This simple observation means that a 
local gauge symmetry cannot strictly speaking be broken spontaneously |29j , and that one should 
only consider manifestly gauge-invariant quantities. 

Nevertheless, it is possible to calculate the effective potential, if one first fixes the gauge. That 
means imposing a constraint, which breaks the gauge invariance, but does not change any gauge- 
invariant quantity. Then, one can calculate the effective potential and use it to study the phase 
transition. This may give non-zero values to non-gauge- invariant quantities such as (</>), but one 
should keep in mind that this is just a non-physical consequence of the gauge fixing. 

Even after the gauge fixing, the calculation of the effective potential is difficult in gauge field 
theories, and the reason is that the gauge bosons are massless in the symmetric phase. If the 
field in the one-loop Feynman diagram that contributes to the effective potential is massless, the 
diagram diverges and the calculation breaks down. This infrared problem can be partly solved by 
using a trick called resummation, but even that only works when the scalar coupling A is weak 
compared to the gauge coupling g 3ft. In that case, there is an extra cubic contribution g 3 Tcj) 3 
to the effective potential. There is, therefore, a range of temperatures in which the potential has 
two minima, one at zero and one at a non-zero value (see Fig. Et>). Furthermore, the values of 
a gauge-invariant quantities such as (|</>| 2 ) are different in these two minima, and therefore they 
actually correspond to physically different states. 

When the minimum that corresponds to non-zero (4>) becomes the global minimum, we have 
exactly the same situation as in Guth's original inflationary model. The symmetric vacuum 
becomes metastable, and bubbles of the true, broken vacuum start to nucleate, in perfect analogy 
with bubbles of vapour in boiling water. A transition like this is known as a first-order phase 
transition. This does not always lead to inflation, though, because in many cases the energy of 
the thermal fluctuations dominates over the potential energy. 

On the other hand, when the scalar coupling is strong relative to the gauge coupling, the 
perturbative calculation of the effective potential is unreliable. Generally, one will then have to 
use numerical lattice Monte Carlo simulations. It is known that for the Abelian Higgs model there 
is a phase transition, and the two phases can be characterized by whether the magnetic field has 
long-range correlations or, in other words, whether there is Meissner effect or not. 




(16) 




(17) 
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These correlations can be seen in the thermal fluctuations of the magnetic field, which are 
actually nothing but thermal radiation, which we will approximate by blackbody radiation. 2 We 
can think of the fluctuations as a superposition of plane waves with different wave numbers k, each 
of which has a Gaussian probability distribution with variance G(k). In the symmetric phase, this 
is simply equal to the energy of the degree of freedom, and because of classical equipartition of 
energy, we have 

G(k) = T. (18) 

In the broken phase, the Meissner effect starts to suppress long-wavelength fluctuations, and we 
find 

G{k) = iT^' (19) 

where to 7 is the inverse of the photon correlation length, which is also known as the penetra- 
tion depth in superconductivity. Near the transition in the broken phase, the value of m 7 is 
approximately given by 

m 2 «e 2 (T 2 -T 2 ). (20) 

In non-Abelian theories such as the electroweak theory, the Georgi-Glashow model or GUTs, 
the transition actually disappears and becomes a smooth crossover when the scalar self-interactions 
dominate over the gauge coupling 31 J. This may sound odd, because one might have thought 
about using the value of (</>) to distinguish between the two phases, but this quantity is not 
gauge-invariant and does not therefore have a physical interpretation. In fact, because the gauge 
symmetry cannot be broken, there is no qualitative difference between the two phases, and in this 
sense the disappearance of the transition should not be surprising. After all, exactly the same 
happens to the water-vapour transition at a high enough pressure. 



8 Global symmetries — the Kibble-Zurek mechanism 

Let us first consider what happens in theories with global symmetries if a symmetry that was 
initially present gets spontaneously broken. The equilibrium state of the system after the tran- 
sition should in principle correspond to a constant, non-zero value of <f> with some small thermal 
fluctuations around it. However, the choice of <f) is not unique, because it could equally well be 
any point on the vacuum manifold. It was pointed out by Tom Kibble in 1976 that this leads to 
defect formation 

If we think of two points that are far enough from each other, there is no reason to expect 
that they would make the same choice. In the cosmological context, an obvious constraint is that 
any two points must make this choice independently of each other if they are so far from each 
other that even a light signal would not have had time to travel from one to the other. However, 
it is equally obvious that this only gives the upper limit for the distance. For the moment, let us 
just say that there is a critical length scale £, and at longer distances the choice of the vacuum is 
uncorrelated. 

In this case, the system cannot reach the true equilibrium state, because that would require 
making the same choice at infinite distances. Instead, the best it can do to minimize the energy 
is to make the field more or less constant at distances less than £. To get a picture of this, we 
can think that it forms domains of radius £ , inside each of which the field is constant and between 
which it is random. 

Of course, there would not really be well-defined domains, because the field would vary 
smoothly from one domain to another. However, even that is not possible everywhere. Imag- 
ine, for instance, that a U(l) symmetry is broken in two dimensions. Fig. shows a possible 

2 Because we are interested in long wavelengths, the fluctuations are described by the classical Rayleigh- Jeans 
spectrum 1181 . This spectrum cannot be valid at very short wavelengths where k>T, because it would lead to 
the famous ultraviolet catastrophe. This was the reason why Planck postulated that the energy of radiation is 
quantized, which eventually led to the development of quantum mechanics. However, because we are considering 
macroscopic effects, the Rayleigh-Jeans spectrum is sufficient. 
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Figure 6: Vortex formation by the Kibble mechanism. Because the system has only a finite 
time available, it can only be ordered at distances less than a critical value £. Therefore, it forms 
domains of these size, inside each of which the scalar field, depicted by the small arrows, is aligned. 
The system tries to interpolate the field between the domains, but it can only do it if the field 
vanishes somewhere between the domains. The vortex appears where the field vanishes. 

configuration of three neighbouring domains, and one can interpolate the field from 1 to 2, from 
2 to 3 and from 3 to 1, but if one does that, one finds that there is a vortex at the point where 
these domains meet. This means that there is roughly one vortex per domain, and the number 
density of vortices formed in the transition is 

n ~ (21) 

The same estimate applies to strings in three dimensions, except that that number density should 
then be interpreted as the number density per cross-sectional area. In other words, the transition 
produces a string network, and the number of strings that cross any two-dimensional plane is given 
by Eq. {2J. 

It is easy to generalize this argument to monopoles, and that yields a number density per unit 
volume of 

n ~ 1/C 3 . (22) 

Because global monopoles have very strong interactions due to their linear interaction potential, 
this expression is only valid immediately after the transition. Later on, the interaction starts 
to pull monopoles and antimonopoles together, and whenever they meet they annihilate. The 
number density of monopoles would therefore drop very quickly from its initial value in Eq. (|22|) . 
To a somewhat lesser extent, the same is also true for global vortices, which have a logarithmic 
interaction potential. 

The predictions in Eqs. Ij21|l and (|22(l are not very useful unless we know how to calculate £. It 
was argued by Zurek in 1985 [2j that it depends mostly on the critical exponents fi and v defined 
in Eqs. I|l(j|) and l|17f) . For simplicity, we assume that the temperature is decreasing linearly, 

T{t) =(l-±-\ T c . (23) 

The quantity tq characterizes the cooling rate, and is known as the quench timescale. When the 
temperature decreases towards T c , the correlation length £ grows as £(t) ~ (I^I/tq) - ^ but at the 
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Figure 7: Snapshot of a liquid crystal after the phase transition. Vortices are the places where 
the dark regions form a cross. (From Digal et al. |H8| ) 

same time the dynamics of the system becomes slower, r(t) ~ (|*|/tq) — p . Eventually, at some time 
i, the relaxation time r becomes equal to the time that is left before reaching T c , i.e., r(i) = \t\. 
After this time, the system cannot adjust to the change of the temperature fast enough, and falls 
out of equilibrium. This freeze-out time depends on tq as i oc —Tq^^. The correlation length 
cannot grow significantly after this time, and therefore we can, to a good approximation, say that 
it freezes to whatever its value was at time i, and that determines the freeze-out scale £ 0, 

£ = m«T^\ (24) 

A similar, but much simpler, picture can be used for first-order transitions that take place 
through bubble nucleation. In that case, the domains would be replaced by bubbles, and the scale 
£ would be the typical bubble radius at the time when the bubbles coalesce. 

The predictions Eq. (|21(1 and (|24|l have been tested experimentally in a variety of systems, 
ranging from superfluid 4 He jS3 B3] and 3 He [S| ESI and li( l uid crystals (see Fig. I7)l 
to nonlinear optical systems |3!5] and convective fluids |I0]. The results agree mostly with the 
predictions, although systematic tests have not been possible in many cases, because it is typically 
quite difficult to vary the quench time scale tq and to determine the critical exponents fi and v. 

Therefore, it is perhaps better to concentrate on how the defects and antidefects are distributed 
in space [5S|. because that is generally insensitive to details like these. Let us, for simplicity, discuss 
a two-dimensional system with vortices and consider the winding number Nw(R) around a circle 
of radius R. On average, this is of course zero, (N\v(R)) = 0, because there are as many vortices as 
antivortices. However, the variance <Tw(R) = (Nw(R) 2 ) is generally non-zero, and its behaviour 
as a function of R reflects the spatial distribution of vortices. 

If we assume that the locations of the vortices are more or less uniformly distributed, the 
number of vortices plus antivortices inside a circle of radius R is A VO rt = mrR 2 , where n is the 
number density. 

If each of the A vort vortices were just as likely to be positive as negative, Nw(R) would be just 
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a result of a random walk of N VOTt steps of unit length. The variance of that is well known to be 

a w (R) = AT vort cx R 2 . (Random) (25) 

However, if the vortices were produced by the Kibble mechanism, we would expect something 
different. The arc of the circle would cross roughly iVd om ~ R/£ correlated domains. Each time 
we move from one domain to another, the phase angle changes by an essentially random amount, 
and around the whole circle, there are Ad om of these changes. Therefore, we have a random walk 
of iVdom steps, and consequently, 

a w {R) oc N dom cx R. (Kibble) (26) 

More generally, we can follow Ref. [38| and characterize the distribution of the vortices by writing 

<j w {R) oc R iv ". (27) 

The above arguments show that a random distribution would give v a — 1/2, whereas the Kibble 
mechanism would predict v„ — 1/4. These predictions do not depend on any quantities that would 
have to be measured or calculated separately such as the critical exponents. 

The exponent v a was measured experimentally for a phase transition in liquid crystals by Digal 
et al. ISHIj who found v a = 0.26 ± 0.11, in agreement with what the Kibble mechanism predicts. 



9 Gauge symmetries — trapping of thermal fluctuations 

If the symmetry broken in a phase transition is a local gauge symmetry, it is not enough to consider 
the behaviour of the scalar field 0T| . The gauge field has a dynamics of its own, and as Eq. J§Jl 
shows in the case of the Abelian Higgs model, the Higgs field does not even try to be aligned if A 
is non-zero. 

Let us start by discussing the Abelian Higgs model. If we imagine having initially a uniform 
magnetic field B = Bo in the system, we know very well from superconductor experiments what 
happens when we cool it into the broken (superconducting) phase jllj : The magnetic field forms 
a lattice of vortices. Because all the flux must go into the vortices, the number density would 
be n = (e/27r)|Bo|. This happens no matter how slowly we cool the system, and therefore the 
formation of these defects is clearly not described by the Kibble mechanism. 

Imagine now that the magnetic field is not uniform, but that there is a magnetic plane wave 
of some wavelength A. In the absence of magnetic monopoles, the magnetic field lines cannot end, 
and the only way in which a field line loop can disappear is therefore by shrinking to a point. 
However, that cannot happen in less time than what it takes for light to travel from one end of the 
loop to the other. This time is simply the wavelength in units of the speed of light, and therefore 
the lifetime of a fluctuation cannot be shorter than the wavelength. It is actually more convenient 
to use the wave number k = 1/A than the wavelength A. If we define the decay rate "f(k) as the 
inverse of the time it takes for the amplitude of a fluctuation with wavenumber k to halve, this 
result gives a causality bound j(k) < k. 

We can also calculate j(k) more precisely. The behaviour of the magnetic field is determined 
by Maxwell's equations, and in vacuum it would simply oscillate with frequency k. However, in 
any medium, such as the hot elementary particle plasma that filled the early universe, there are 
electric currents. If they obey Ohm's law with conductivity a, long-wavelength modes (k <C cr) 
have a decay rate 

7 (fc) = (28) 

which is much slower than the causality bound. 

If we start with the magnetic plane wave and cool the system into the superconducting phase 
very rapidly, the wave does not have time to decay. Because there can be no magnetic field in 
the broken phase, it forms vortices instead, as illustrated in Fig. [5] The vortices point in the 
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a) 




Figure 8: Vortex formation by flux trapping in the Abelian Higgs model, a) In the symmetric 
phase, there are Gaussian thermal fluctuations of all wavelengths, b) During a rapid quench, 
longest wavelengths A >l/k c are too slow to react and freeze out. Shorter wavelengths disappear, 
leading to a smooth magnetic field, which varies over the distance scale l/k c . c) The broken phase 
is a superconductor, and the frozen-out magnetic field becomes trapped into vortices. 
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same direction as the initial magnetic field, and their number density is higher where the field 
was stronger. In fact, if we looked at the system with a low enough resolution so that we could 
not see individual vortices but could still measure the average magnetic field arising from the flux 
quanta carried by the vortices, we would still see the initial magnetic field. On the other hand, if 
the cooling is slow enough, then the initial wave decays and no vortices are formed. 

So, we have seen that a non-zero initial magnetic field leads to vortex formation, but where 
could that initial magnetic field come from? The answer is thermal fluctuations. We already 
saw that arbitrarily long wavelengths are present in the Rayleigh- Jeans spectrum Ijl8|l of thermal 
radiation. Let us consider a single mode with wave number k, i.e. a single plane-wave component 
of the thermal spectrum. If the corresponding decay rate 7(A) is slower than what it would have 
to be in order for the mode to stay in equilibrium, the mode freezes out and keeps its initial 
equilibrium amplitude, forming vortices as we saw above. 

If we know how the temperature changes, we know from Eq. I|19|) . what the equilibrium ampli- 
tude of any given mode is and therefore how quickly the mode would have to to decay in order to 
adjust to this change. It is obvious that there are always going to be some very long wavelengths 
that are unable to do it. Simply from causality, we obtain the lower bound 

e 2 T 2xl/3 



and if Ohm's law is valid, we have 



fc c >(-^) , (29) 



e 2 T 2 X l/4 



a . (30) 
Q J 

Now, we will have to find out what this means in terms of the number density of vortices. By 
definition, the modes with k < k c freeze out, and those with k > k c equilibrate. This magnetic 
field is squeezed into vortices, but let us ignore that for a moment, and assume that all modes 
with k > k c simply disappear. This means that the magnetic field configuration is a superposition 
of waves with wavelength longer than l/fc c , and if we look at the system at distances shorter 
than that, we simply see a uniform magnetic field. Eventually, all this flux turns into vortices (see 
Fig. ISJ;). If there was initially a flux $(-/?) through a circle of radius R< l/k c , it will form <I>(i?) / <&o 
vortices, all with the same sign. The number density n can then be calculated by dividing this by 
the area irR 2 of the circle. 

On the other hand, if R>l/k c , the Fourier modes with k > k c only affect the spatial distri- 
bution of the flux inside the curve, but do not contribute to the net flux. Consequently, the final 
flux Q(R) after the transition is the same as it was in the symmetric phase before the transition. 
This can be calculated from G(k) or estimated using physical arguments, and the result is |41| 

$(i?) w VTR. (31) 

If we now choose R = l/k c , we are on the borderline of these limiting cases. We can say 
that the flux is given by Eq. (|31|) . and that the whole flux is turned into vortices of equal sign. 
Therefore, we know what <i>(l/fc c ) is, and we can calculate the number density of vortex lines, 



$(i?)/$ 

n ~ 



ttR 2 

Using Eq. (|30|) . we can write this as 



R=l/k a 



V^Tfcf. (32) 



7/4 T 5/4( M 3/S . (33) 



n ss e 



The two-dimensional analogue of this result has been tested in numerical simulations 021 > an d 
seems to work well (see Fig. 0). There have also been attempts to test the result experimentally 
with superconductors (see Fig. HT1|l [4"3l EH 0Sj. However, the actual experiments are usually 
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Figure 9: A snapshot from a numerical simulation of the Abelian Higgs model. It shows the 
magnetic field configuration after a quench into the broken phase. Black and white correspond to 
magnetic fields pointing up and down, respectively. The clustering of vortices is obvious. (From 
Stephens et al. |32"]1. 

done with two-dimensional films. This is not equivalent to the fully two-dimensional setup used 
in the simulations, because in the experimental setup, the magnetic field extends outside the 
film. Therefore, the magnetic field behaves as in three dimensions and the scalar field as in two 
dimensions. This situation was discussed briefly in Ref. |41|. and a more detailed analysis is given 
in Ref. gBJ. 

However, just like the Kibble mechanism, this scenario can also give predictions that are 
independent of dynamical parameters such as the conductivity. As we saw, the flux <&(-/?) through 
a circle with R = l/k c turns into vortices of equal sign, and therefore, if it is greater than 
<J>0i we should see clusters of equal-sign vortices. Based on Eq. i|32|) . we would expect roughly 
N cl « {e 2 Tlk c f/ 2 vortices per cluster. Because we can make k c arbitrarily small by cooling the 
system slower, we should be able to make these clusters if we can control the cooling rate. 

If N c \ is large enough, it is much more likely that a randomly picked vortex is inside a cluster 
than near it boundary. Therefore, if we calculate the quantity ow{R) defined in Section |H1 all the 
vortices contributing to Nw(R) are of the same sign as long as R is less than the cluster size l/k c . 
Therefore Nw (R) is proportional to its area, and consequently, 

aw(R) oc i? 4 , (Flux trapping) (34) 

at short distances. In other words, v a = 1. 

Furthermore, if we can measure both the number density n and the cluster size N c i, the 
combination nN^ should be independent of k c and can therefore be used to test the scenario. On 
the other hand, if one wants to learn more about the dynamics of the phase transition, one would 
measure that combination n/N c \, which would give the freeze-out scale k c . 

The result in Eq. l|3"2"|) is proportional to the square root of temperature and would therefore 
vanish as T goes to zero. This is important, because the transition that ended hybrid inflation 
would have taken place at zero temperature. Would defect formation then have happened accord- 
ing to the Kibble mechanism, or would the gauge fields still have played some role? In the above 
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Figure 10: The distribution of magnetic flux in an array of superconducting rings after a rapid 
phase transition. The similarity with Fig. [5] is striking, and the authors report that their re- 
sults show that thermal fluctuations play an important role. The result is still not a conclusive 
proof of the theory discussed in Section [5] because a ring array is in many ways different from a 
superconducting film or indeed a three-dimensional space. (From Kirtley et al. |44p 

discussion, we have only considered classical thermal fluctuations, and it might be possible that 
quantum fluctuations play a similar role at zero temperature. More work is needed to find out 
whether that is really the case. 

I have argued recently in Ref. 07J that the formation of magnetic monopoles in the GUT 
phase transition may take place in a somewhat similar way. There are several extra complications, 
though. In particular, as was already mentioned, there may actually be no true phase transition 
at all, just a smooth crossover from one phase to another. One may then ask whether monopoles 
are formed at all. 

It is, however, easy to see that the theory cannot stay in equilibrium if it is cooled at a constant 
rate. There are no long-range correlations in the symmetric phase of the theory, and it is generally 
known that any correlation length has to be of the order of (<? 2 T) _1 or shorter ^Sji where g is again 
the gauge coupling constant. On the other hand, at a very low temperature we should recover 
normal electrodynamics, which has a massless photon. Correspondingly, the magnetic field has an 
infinite correlation length. 

Admittedly, it is not obvious what "magnetic field" means in the symmetric phase, but no 
matter how we define it, we have a correlation length £b that starts from a microscopic value 
(g 2 T) _1 and grows to an infinite value. Analogously with Kibble's argument, we can say that 
£b cannot grow arbitrarily fast, because causality alone restricts its growth rate to be less than 
the speed of light. It is therefore inevitable that if we keep on cooling the system, it falls out of 
equilibrium at some point and the magnetic correlation length freezes to a finite value, which we 
denote by £b- 
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The next key observation is that the finiteness of the correlation length is due to screening. 
This is analogous to the Debye screening of electric field by electric charges 05] ; an d therefore it 
implies a presence of magnetic charges. Thus, because the correlation length is necessarily finite 
due to causality, magnetic monopoles must have been formed. 

These arguments can be made more quantitative by noting that one can define a magnetic 
charge that is conserved even in the symmetric phase. This charge behaves in a similar way to 
the magnetic field in the Abelian Higgs model. In the symmetric phase, there are long-wavelength 
thermal fluctuations, which should disappear in the broken phase, but again, causality implies 
that fluctuations whose wavelength exceeds some critical value do not have time to decay. They 
freeze out and form vortices. This critical wavelength turns out to be exactly the value £b to 
which the correlation length of the magnetic field freezes. One finds that the number density of 
monopoles formed at the transition would be roughly 07] 



As for vortices, one also find that if £b is long enough, clusters of y <? 2 T£b monopoles are formed. 



So far, we have been discussing continuous (or smooth) phase transitions, but discontinuous 
first-order transitions are possible both in Abelian and non-Abelian theories, and in general they 
take place if the Higgs self coupling constant A is less than the square of the gauge coupling 
constant. 

In a first-order transition, bubbles of the broken phase are nucleated, and they grow and 
coalesce. We can safely assume that there are no defects inside the bubbles. Because magnetic 
field and charge are conserved in Abelian and non-Abelian cases, respectively, they will be trapped 
in the gaps between the bubbles. Therefore, we would expect that monopoles are formed in well- 
localized clusters. In the Abelian case, vortices with multiple windings are stable if there is a 
first-order transition, and therefore we would expect thick and massive vortex lines, which carry 
much higher fluxes than in the continuous case |5l)j . 

10 Defects in cosmology 

There is no conclusive evidence for the existence of topological defects in the universe, but there 
is strong evidence against certain types of topological defects. First of all, domain walls are ruled 
out unless they, for some reason, have an extremely low energy density 0]. The reason for this 
is that the mass of a domain wall is proportional to its area and is therefore enormous for a wall 
that extends through the whole universe. 

Global strings or monopoles could arise from a spontaneous breakdown of a continuous global 
symmetry, but as we saw in Section |3] that would predict the existence of massless Goldstone 
bosons. No such particles have been observed, but they may still exist if they only interact 
extremely weakly with other forms of matter. 

Otherwise, we are left to consider gauged defects. Cosmic strings have been studied intensively, 
because it was thought they could explain the primordial density perturbations that gave rise to 
the large-scale structure of the universe and the temperature fluctuations of the cosmic microwave 
background (CMB) radiation. However, recent CMB measurements ^H] show that this is not the 
case, and that these perturbations probably originated in quantum fluctuations during inflation. 
This observation does not by any means rule out the existence of cosmic strings, and people 
are still looking for evidence for or against them in the CMB radiation and in other astronomical 
observations. One way of finding strings is by gravitational lensing. Because strings are so massive, 
they would bend the light coming from galaxies behind them and therefore we would see two copies 
of the same galaxy. A possible such observation of a cosmic string has been reported recently |51| . 

In principle, magnetic monopoles would be a natural candidate for the missing dark matter, 
which is known to make up about 85% of all the matter in the universe [TBI ■ The y have been 




(35) 
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searched for in experiments, but with no success [52] ■ This seems to rule out the possibility of 
superheavy magnetic monopoles with GUT scale masses being the main component of the dark 
matter. There is also a more severe problem: The number of magnetic monopoles produced at 
the GUT phase transition has been estimated to be so high that their mass would have caused 
the universe to collapse very soon after the transition |531 154] . This is known as the monopole 
problem. 

In the standard Big Bang scenario with no inflation, we would assume that the temperature 
was initially around the Planck scale 10 19 GeV. When the universe was expanding and cooling 
down, it went though the GUT transition at around 10 16 GeV. Because all GUTs predict the 
existence of magnetic monopoles and they would have been formed at this transition, this scenario 
suffers from the monopole problem. 

The monopole problem is solved by inflation if the reheat temperature was low enough so that 
the GUT symmetry was never restored. So far, there are very few other observational constraints 
on the reheat temperature, and more or less the only other constraint comes from the matter- 
antimatter asymmetry. There is overwhelming observational evidence that there is practically no 
antimatter but a significant density of matter in the universe. It seems impossible to generate 
this asymmetry after the electroweak phase transition [HE], which took place at around 100 GeV, 
and therefore the reheat temperature should be high enough for the electroweak symmetry to be 
restored. This gives the rough bounds, 

100 GeV < T rh < 10 16 GeV. (36) 

If there were other phase transitions, cosmic strings may have been produced provided that the 
reheat temperature was above the corresponding critical temperature. 

There could, however, be extra complications. First of all, it is not at all clear, how well the 
scenarios that assume an equilibrium initial state can really describe phase transitions that take 
place in a rapidly expanding universe only a short time after the end of inflation. The assumption 
of thermal equilibrium seems particularly questionable at long distances, between points that have 
only barely come within each other's sight. 

In many inflationary models, reheating can be a very violent non-equilibrium process, in which 
case it is often called preheating jSH]- It may then be possible to produce topological defects |57l 
1581 I59|R)U| even if the reheat temperature is well below the critical temperature. In some cases, it 
is possible to describe this as a "non-thermal" symmetry restoration. It is not known how well the 
scenarios discussed in Sections [S] and |5] can describe defect formation from preheating. In any case, 
it seems that we may have to rethink the upper bound in Eq. H3fi[) ■ The same applies to the lower 
bound, because similar effects may also create the matter-antimatter asymmetry [21 EH EH! ■ 

There is also a third possibility in hybrid inflation l)15[l. Because inflation ends at a breakdown 
of a symmetry, there is always a phase transition at the end of inflation, and that would typically 
lead to defect formation If 

we want to avoid domain walls and Goldstone bosons, 
the waterfall field x wm have to be a multi-component field with a gauge symmetry. Because 
the universe is at zero temperature at the time of the phase transition, it is not yet clear if the 
scenario discussed in Sectional is applicable. 

11 Conclusions 

Theoretically, it seems plausible that topological defects should have been formed in the early 
universe. If signs of them are observed, that will give an enormous boost to our understanding 
of the universe, but the failure to find them can also tell us many things. Most significantly, 
the absence of the high number of monopoles predicted by the theory led to the development 
of the theory of inflation, which has become the standard paradigm in cosmology because of 
other supporting observations. The apparent absence of defects can also help us constrain certain 
parameters of the inflationary models, which would otherwise be undetermined. This will become 
more important in the future, when we start to home in on a particular class of inflationary 
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models using data from other observations. It also has significance in particle physics, because it 
constrains possible high-energy extensions of the Standard Model of particle physics. 

For all this it is important to have a good understanding of how defects are formed. Currently, 
there are theories that describe this process in transitions that start from thermal equilibrium 
in either global or gauge field systems. In the global case, some of the predictions have been 
confirmed in condensed matter experiments, and it seems that the gauge theory predictions will 
be tested with superconductors in the future. 

However, the transitions in the early universe may have taken place far away from thermal 
equilibrium, and much more work is still needed before we can really claim that we understand 
how possible cosmological defects would have been formed. 
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